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ABSTRACT 

We investigate the relationship between the mass of a globular cluster core and the 
sizes of its various stellar populations in a sample of 56 globular clusters. The number 
of core red giant branch stars is found to scale sub-linearly with core mass at the 
3-(T confidence level, whereas the relation is linear to within one standard deviation 
for main-sequence turn-off and sub-giant branch stars. We interpret our results as 
evidence for a surplus of red giant branch stars in the least massive cluster cores 
which is not seen for main-sequence turn-off and sub-giant branch stars. We explore 
various possibilities for the source of this discrepancy, discussing our results primarily 
in terms of the interplay between the cluster dynamics and stellar evolution. 

Key words: globular clusters: general - stars: statistics - stellar dynamics - stars: 
evolution. 



1 INTRODUCTION 

Studying the radial distributions of the various stellar popu- 
lations (red giant branch, horizontal branch, main-sequence, 
etc.) found in globular clusters (GCs) can provide useful 
hints regarding their dynamical histories. As clusters evolve, 
they are expected to undergo relatively rapid mass stratifi- 
cation as a consequence of two-body relaxation, with the 
heaviest stars quickly sinking to the central cluster regions 
(|Spitzeilll987h . The shorter the relaxation time (typically 
evaluated at the half-mass radius), the quicker this process 
occurs. Clusters tend towards dissolution as two-body relax- 
ation progresses and they lose mass due to stellar evolution 
and the preferential escape of low-mass stars. External ef- 
fects like tidal perturbations, encounters with giant molecu- 
lar clouds, and passages throug h the Galactic disk serye only 
to speed up the pr ocess (e.g. iBaumgardt fc Makinol l2003l : 
iKiipper et al.ll2008l ). Stellar evolution complicates this oth- 
erwise simple picture of GC evolution, however. Stars are 
expected to change in size and lose mass as they evolve, 
often dramatically, and this could significantly impact the 
outcomes of future dynamical interactions with other stars. 
For instance, a typical star in a GC is expected to expand by 
up to a few orders of magnitude as it ascends the red giant 
branch (RGB) and will shed up to a quarter of its mass upon 
evolving fr om the tip of the RGB to the horizontal-branch 
(HB) (e.g. lCaloi fc D'Antona.200a : .Lee et al...l994 ). 
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Red giant branch stars have been reported to be defi- 
cient i n the cores of some Milky Way (MW) GCs. For in- 
stance, iBailvnl (Il994l ) found that the morphology of the gi- 
ant branch in the dense core of 47 Tuc differs markedly from 
that in the cluster outskirts. In particular, there appear to 
be fewer bright RGB stars in the core as well as an enhanced 
asymptotic giant branch (AGB) sequence. While a similar 
deficiency of bright giants has been observed in the cores of 
the massive GCs NGC 2808 and NGC 2419, better agree- 
ment between the observations and theoretical luminosity 
functions obta ined with the Victor ia-Rcgin a isochrones was 
found for M5 JSandguist fc Martel,.2007. : Sandauist fc Hesj 
I2OO8I) . ISandauist fc Martell (|2007l ) speculate that the giant 
star observations in NGC 2808 could be linked to its unusu- 
ally blue horizontal branch if a fraction of the stars near the 
tip of the RGB experience sufficiently enha nced mass loss 
that t hey leave the RGB early. Alternatively. iBeer fc DaviesI 
(|2004l ) suggest that RGB stars could be depleted in dense 
stellar environments as a result of collisions between red gi- 
ants and binaries. 

While stellar populations have been studied and com- 
pared on an individual cluster basis, a statistical analysis in 
which their core populations are compared over a large sam- 
ple of GCs is ideal for isolating trends in their differences. 
Though a h andful of studies of this nature have been per- 
formed (e.g. IPiotto et al.ll2004h . we present an alternative 
method by which quantitative constraints can be found for 
the relative sizes of different stellar populations. Specifically, 
a cluster-to-cluster variation in the central stellar mass func- 
tion can be looked for by comparing the core masses to the 
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sizes of their various stellar populations. Since stellar evolu- 
tion is the principal factor affecting their relative numbers in 
the core, we expect the size of each stellar population to scale 
linearly with the core mass. If not, this could be evidence 
that other factors, such as stellar dynamics, are playing an 
important role. In this paper, we present a comparison of the 
core RGB, main-sequence (MS) and HB populations of 56 
GCs. In particular, we use star counts for each stellar pop- 
ulation to show that RGB stars are either over-abundant 
in the least massive cores or under-abundant in the most 
massive cores, and that this effect is not seen for MS stars. 
We present the data in Section [5] and our methodology and 
results in Section |3] In Section [l] we discuss the implica- 
tions of our results and explore various possibilities for the 
source of the observed discrepancy between RGB and MS 
stars. Concluding remarks are presented in Section [5] 



2 THE DATA 

Colour-magn i tude diagrams (CMDs) taken from 
IPiotto et al] (|2002l ) are used to obtain star counts for 
the RGB, HB, MS and blue straggler (BS) populations in 
the cores of 56 GCs. We apply the same selection criterion 
as outlined in Leigh, Sills & Knigge (2007) to derive our 
sample as well as to define the location of the main-sequence 
turn-off (MSTO) in the (F439W-F555W)-F555W plane. 
An example of this selection criterion, applied to the CMD 
of NGC 362, is s hown in Figure [1] We include all s t ars in 
the lPiotto et al.l (|2002l 'l database. Since iPiotto et al.l (|2002l ') 
took their HST snapshots with the centre of the PC chip 
aligned with the cluster centre, a portion of the cluster core 
was not sampled for most GCs. We have therefore applied 
a geometric correction factor to the star counts in these 
clusters in order to obt ain numbers that arc representative 
of the entire core (|Leigh. Sillsfc Kniggel |2007| , |2008j). The 
total number of stars in the core is found by summing over 
all stars brighter than 1 mag below the MSTO and then 
multiplying by the appropriate geometric correction factor. 

Errors on the number of stars for each stellar population 
were calculated using Poisson statistics. Core radii, distance 
moduli, extinction corrections, central luminosity densities 
and central surface brightnesses were ta ken from the Harris 
Milky Way Globular Cluster catalogue (jHarrisI [l996h . Cal- 
ibrated apparent magnitudes in the F555W, F439W and 
Johnson V bands were taken from iPiotto et al] (|2002l ). 



3 RESULTS 

This paper focuses on the core RGB, MS and HB popula- 
tions of 56 GCs, comparing their numbers to the core masses. 
Note that we are focusing on the total number of stars in 
the core as a proxy for the core mass instead of the total 
luminosity in the core in order to avoid concerns regarding 
cluster-to-cluster variations in the central stellar mass func- 
tion and selection effects. Given that a single bright HB star 
can be as luminous as 100 regular MS stars, a small surplus 
of bright stars could have a dramatic impact on the total 
luminosity. Therefore, the total number of stars in the core 
is a more direct and reliable estimate for the core mass than 
is the core luminosity. 
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Figure 1. Colour-magnitude diagram for NGC 362 in the 
(F439W-F555W)-F555W plane. Boundaries enclosing the se- 
lected RGB, HB and MSTO populations are shown. 

Upon plotting the logarithm of the number of core RGB 
stars versus the logarithm of the total number of stars in the 
core and performing a weighted least-squares fit, we find a 
relation of the form: 

log{NRGB) = (0.89 ± 0.03) log(Ar,<,,,/10^) + (2.04 ± 0.02) 

(1) 

The sub-linear slope is either indicative of a surplus of RGB 
stars in the least massive cluster cores or a deficiency in the 
most massive cores. Errors for lines of best fit were found 
using a bootstrap methodology in which we generated 1,000 
fake data sets by randomly sampling (with replacement) 
RGB counts from the observations. We obtained lines of 
best fit for each fake data set, fit a Gaussian to the sub- 
sequent distribution and extracted its standard deviation. 
In order to avoid the additional uncertainty introduced into 
our RGB number counts from trying to distinguish AGB 
stars from RGB stars, as well as the difficulty in creating a 
selection criterion that is consistent from cluster-to-cluster 
when including the brightest portion of the RGB, stars that 
satisfy the RGB selection criterion shown in Figure [1] are re- 
ferred to as RGB stars throughout this paper. Note that it 
is the brightest portion of the RGB that should be the most 
affected by stellar evolution effects such as mass-loss. If we 
extend our selection criterion to include the entire RGB, 
however, our results remain unchanged. 

Interestingly, MS plus sub-giant branch stars (hereafter 
collectively referred to as MSTO stars, the selection criterion 
for which is shown in Figure 1) show a more linear relation- 
ship than do RGB stars and appear to dominate the central 
star counts. If we count only those stars having a F555W 
mag within half a magnitude above and below the turn-off, 
we obtain a relation of the form: 

log(iVMSTo) = (1.02 ± 0.01) log(iV,„,e/10^) + (2.66 ± 0.01) 

(2) 
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Figure 2. The number of RGB (red circles), MSTO (black 
squares) and HB (blue triangles) stars found within the cluster 
core plotted versus the total number of stars in the core brighter 
than 1 mag below the MSTO, along with the corresponding lines 
of best fit for the RGB and MSTO samples. 



A nearly identical fit is found when counting only those stars 
having a F555W mag between the turn-off and one magni- 
tude fainter than the turn-off. 

We also tried plotting the logarithm of the number of 
core helium-burning stars (labeled HB in Figure [1]) versus 
the logarithm of the number of stars in the core, yielding a 
relation of the form: 

log(iVHB) = (0.91±0.10) log(iV,„,,/10^)-f (1.58±0.05) (3) 

Note the large uncertainty associated with the fit, indicating 
that the slope is consistent with both those of the RGB 
and MSTO samples. We will discuss this stellar population 
further in Section 

The number of MSTO, RGB and HB stars are shown 
as a function of the total number of stars in the core in 
Figure (2] Interestingly, the blue stragglers in our sample 
also scale sub-linearly with core mass, albeit more dramat- 
ically, obeying a relation o f the form Nbs ^ M°o'^g'^'' °* 
l|Knigge. Leigh fc Sill3l2009l ). Note that Ncore can be used 
interchangeably with Mcore- In this case, we obtain a fit of: 



log(iVss) = (0.47±0.06)log(iV,o,e/10^) + (1.22±0.02) (4) 

In an effort to explore the infiuence of selection effects, 
we re-did our plots having removed from our sample clus- 
ters denser than log p > 10^ Lq pc~^ since we are the most 
likely to be under-counting stars in the most crowded clus- 
ter cores where blending of the stellar light is the most se- 
vere. This cut also removes from our sample the post-core 
collapse (PCC) clusters for which the core radii are poorly 
defined since King models are known to provide a poor fit 
to the observed surface brightness profiles in these clusters. 
Similarly, we applied a cut in the central surface brightness, 
removing from our sample clusters satisfying Eo < 15.1 V 



mag arcsec"^. Finally, since clusters having both high sur- 
face brightnesses and small cores are the most likely to suffer 
from selection effects, we also tried adding to the aforemen- 
tioned cut in Eo a cut in the angular core radius, removing 
clusters with Tc ^ 0.05'. In all cases, the sub-linear power-law 
index reported for the RGB remains unchanged to within 
one standard deviation of our original result. Selection ef- 
fects do not appear to be the source of the observed sub- 
linearity, though it is clear that its effects must properly be 
accounted for in future studies. 

In order to assess the effects of age-related cluster-to- 
cluster variations in the stellar mass function, as well as to 
test our assumption that the number of stars in the core 
provides a reliable estimate for the core mass, we have ob- 
tained MSTO masses for most of the G Cs in our sample . 
We fit theoretical isochrones provided in I Pols et ahl (| 19981 ) 
to the cluster CMDs, using the bluest point along the MS of 
a given isochrone as a proxy for the MS TO mass. Isochrones 
were calculated using the metalli cities of Piotto et al.l (|2002| ) 
and cluster ages were taken from lDe Angeli et al.l ( 20051 ) us- 
ing the Zinn & West (1984) metallicity scale. Core masses 
were estimated by multiplying the mass corresponding to 
the MSTO {riiMSTo) by the number of stars in the core 
brighter than 1 mag below the turn-off. This is a reasonable 
assumptio n given the very small dispersion in the ages of 
MW GCs l|De Angeh et al.ll2005l ) and the fact that we are 
only considering stars brighter than 1 mag below the TO. 
Consequently, the range of stellar masses upon which we are 
basing our number counts is very small. Our results remain 
entirely unchanged upon using Mcore ^ i^core'mMSTO as a 
proxy for the core mass instead of pure number counts. 

In order to further check the sensitivity of our results to 
our estimate for the core masses, we re-did all plots shown 
in Figure [2] using various approximations for the total core 
luminosity instead of pure number counts. Core luminosi- 
ties are calculated in the Johnson V band directly from the 
stellar fluxes which are summed over all stars in the core 
and then multiplied by the appropriate geometric correc- 
tion factor. We also adopted Lcore ~ ^Tvrlpo, wher e po is 
the ce ntral luminosity density in Lq pc~'^ taken from lHarrij 
(1 19961 ). Additionally, since the number of core RGB stars is 
in reality a projected quantity, we tried plotting N_rgs ver- 
sus Lcore = Trr^Eo, where Eo is the central surface bright- 
ness in Lq pc~^, so that we are consistently comparing two 
projected quantities. Finally, we can adopt slightly more re- 
alistic estimates for the total core luminosity by integrating 
over King density profiles. We fit s ingle-mass King mod- 
els ca lculated using the method of ISigurdsson fc PhinnevI 
(| 19951 ) to the surface brightness profiles of the majority of 
the cl usters in our sample using the concen tration parame- 
ters of iMcLaughhn fc va n den Marell (|2005h and the central 
luminosity densities of iHarris ( 199^yrWe then integrated the 
derived luminosity density profiles numerically in order to 
estimate the total stellar light contained within the c ore. Af- 
ter rem oving clusters labelled as post-core collapse in lHarrij 
( jl996l ) for which King models are known to provide a poor 
fit, we once again compared the integrated core luminosi- 
ties to the number of RGB stars in the core. For all four of 
these estimates for the total core luminosity, we find that our 
fundamental results remain unchanged, with the power-law 
index for RGB stars remaining sub-linear at the 3-a level. 
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Therefore, we conclude that our result is robust to changes 
in choices of cluster and population parameters. 



4 DISCUSSION 

We have shown that the number of RGB stars in globular 
cluster cores does not directly trace the total stellar popula- 
tion in those cores. In particular, the number of RGB (but 
not MSTO) stars scales sub-linearly with core mass at the 3- 
o" level. Given that the MS lifetime is expected t o be a facto r 
of 10-100 longer than that of the RGB sample ijlbenlligQlT l. 
the ratio '^^msto /'^rgb indicates that the relative sizes of 
these stellar populations are in better agreement with the 
expectations of stellar evolution theory in the most massive 
cores. This suggests that our results are consistent with a 
surplus of RGB stars in the least massive cores. We discuss 
below some of the key considerations in understanding the 
evolution of GC cores and the stars that populate them in 
an effort to explain our result. 

4.1 Stellar evolution 

Could this trend be a reflection of a stellar evolution pro- 
cess? The evolution and distribution of stellar populations 
can be thought of as the sum of many single stellar evolution 
tracks, which depend only on a star's mass and composition. 
Since there is no relatio n between a cluster's mass and its 
metallicity (iHarrisI [19961) and the dispersion in the r elative 
ages of MW GCs is quite small (|De Angeli et al.ll2005l ). there 
is no reason to expect the RGB lifetime to depend on the 
cluster mass. On the other hand, recent studies suggest that 
the chemical self-enrichment of GCs during their early evolu- 
tionary stages could help to explain man y of the population 
differ ences observed among them (e.g. ICaloi fc D'Antonal 
I2OO7I ). In particular, many of the most massive GCs are 
thought to be enriched in helium and this is expected to re- 
duce the time scale for stellar evolution (e.g. I Romano et al.l 
I2OO7I ). While this scenario predicts a deficiency of RGB stars 
in the most massive cores, it would contribute to depressing 
the slope of the RGB sample relative to that of the MSTO 
population in Figure (2] 

4.2 Single star dynamics 

Two-body relaxation is the principal driving force behind 
the dynamical evolution of present-day GCs, slowly steer- 
ing them towards a state of increased mass stratification 
as predominantly massive stars fall into the core and typi- 
cally low-mass stars are ejected via dynamical enc ounters. 
The relaxation time increases with the cluster mass l|Spitzeil 
Il98it ) and the variance in the relative ages (and hence MSTO 
meisses) of MW GCs is quite small yZ)c Angeli et al. 2005). 
Therefore, it is the least massive clusters that should show 
signs of being the most dynamically evolved. This assumes 
that cluster-to-cluster variations in the initial mass function 
and the degree of initial mass segregation are small. In gen- 
eral, however, proportionately fewer massive stars should 
have had sufficient time to migrate into the most massive 
cores, while fewer low-mass stars should have been ejected 
out. While qualitatively correct, this effect should contribute 
little to the observed difference between the core RGB and 



MSTO popu lations since RGB stars are only slightly more 
massive (e.g. iDe Marchi fc PulondlioOTl ). 

Stars expand considerably as they ascend the RGB. 
Both the increase in coUisional cross-section and the change 
in the average stellar density could have an important bear- 
ing on the outco mes of dynam ical interactions involving 
RGB stars. Indeed. iBaihml |l994l ') suggests that interactions 
between giants and other cluster members in the core could 
strip the outer envelope of the giant before it has a chance 
to fully ascend the RGB. Since our adopted RGB selection 
criterion does not include the brightest giants, we are only 
considering giants that are larger than MSTO stars by a fac- 
tor of ~ 10 (|lbenlll99ll ). This small degree of expansion will 
have only a minor effect on the collision rate. Any scenario 
that relies on dynamical encounters to explain a depletion 
of RGB stars should be operating in very dense cores. Our 
results are consistent with some of the densest clusters in 
our sample having a surplus of giants, however. 

4.3 Binary effects 

Stripping of the envelopes of large stars could also be me- 
diated by a binary companion a s the expanding giant over- 
fills its Roche lobe (|Bailvnlll99j ). While this process should 
preferentially occur in the centres of clusters where bina- 
ries will congregate as a result of mass segregation, two- 
body relaxation progresses more slowly in the most mas- 
sive clusters. Binaries should therefore sink to the cluster 
core more quickly in the least massive clusters, contribut- 
ing to an increase in the core binary fraction at a rate that 
decreases with increasing cluster mass. Observational evi- 
dence has been fou nd in support of this, m ost notably by 
ISollima et al] (120071 ) and lMilone et al.1 l|2008l ) who found an 
anti-correlation between the cluster mass and the core bi- 
nary fraction. Any mechanism for RGB depletion that relies 
on binary stars should therefore operate more efficiently in 
the least massive cores where the binary fraction is expected 
to be the highest. Our results are consistent with a surplus of 
giants in the least massive cores, however. This therefore ar- 
gues against a binary mass-transfer origin for RGB depletion 
in massive GC cores. For similar reasons, it seems unlikely 
that our result can be explained by collisions between RGB 
stars and binaries. If, on the other hand, RGB stars are more 
commonly found in binaries than are MS stars, perhaps as a 
result of their larger cross-sections for tidal capture, binary 
stars could still be contributing to the observed trend. Note 
that in the cluster outskirts where the velocity dispersion 
has dropped considerably from its central value, individual 
encounters are more likely to result in tidal capture. Since 
mass segregation should deliver binaries to the core faster 
in the least massive clusters, a larger fraction of their RGB 
stars could have hitched a ride to the core as a binary com- 
panion. However, both the average half-mass relaxation time 
of MW GCs and the RGB lifetim e tend to be on the order 
of a Gyr (|Harrislll996l : llbenlll99lh . This does not leave much 
time for giants to be captured into binaries and subsequently 
fall into the core before evolving away from the RGB. 

4.4 Core helium-burning stars 

The fit for the HB sample is consistent with those of the 
RGB and MSTO samples at the 3-sigma level so that we 
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are unable to draw any reliable conclusions for this stellar 
population. The high uncertainty stems from a number of 
outlying clusters. Selection effects and contamination from 
the Galaxy are likely to be playing a role in this, in ad- 
dition to our formulaic selection criterion which may not 
be as suitable to the varying morphology of the HB as it 
is to other stellar populations. That is, the creation of a 
purely photometry-based cluster-independent selection cri- 
terion may not be possible for HB stars. Given that stellar 
evolution effects are expected to be the most dramatic at 
the end of the RGB lifetime, an interplay with the cluster 
dynamics could also be contributing. In particular, if the 
central relaxation time is shorter than the HB lifetime, sig- 
nificant numbers of HB stars could be ejected from the core 
via dynamical encounters as a result of having lost around a 
quarter of their mass upon evolving off the tip of the RGB. 
Moreover, since stars expand considerably as they ascend 
the RGB, many of the dynamical arguments presented in 
Section 14.21 may more strongly affect the size of the HB 
sample if they are the direct evolutionary descendants of 
RGB stars. Since at most a handful of studies have been 
performed c omparing the radial H B and RGB distributions 
in GCs (e.g. Ilannicola et al.ll2009l '). more data is needed be- 
fore any firm constraints can be placed on the source of the 
poor fit found for the HB sample. 

4.5 An evolutionary link with blue stragglers? 

The addition of a small number of extra RGB stars to ev- 
ery cluster is one way to account for the observed sub-linear 
dependence on core mass since the fractional increase in the 
size of the RGB population will be substantially larger in the 
least massive cores. In log-log space, the result is a reduction 
of the slope. Since bl ue stragglers wi ll evolve into RGB and 
eventually HB stars (jSills et al.ll2009l ). evolved BSs could be 
the cause of a surplus of RGB (and possibly core helium- 
burning) stars in these clusters. This scenario also predicts 
that MSTO stars should scale slightly more linearly with 
core mass since there should be a smaller contribution from 
evolved BSs, as we have shown. Given the fits for the RGB 
and BS samples presented in Section|3]and their correspond- 
ing uncertainties, we find that the addition of evolved BSs to 
the RGB populations could inflate the slope enough that the 
dependence on core mass becomes linear. Upon subtracting 
the BS sample from the RGB sample, we find that the new 
fit is consistent with being linear: 

log(iVHGs) = (0.94 ± 0.04) log(iVeore/10^) + (1.97 ± 0.02) 

(5) 

The slope becomes larger if we have under-estimated the 
number of BSs, perhaps as a result of selection effects, our 
adopted selection criterion or a larger population size having 
existed in the past. 



5 SUMMARY 

In this paper, we have performed a cluster-to-cluster com- 
parison between the number of core RGB, MSTO & HB 
stars and the total core mass. We have introduced a tech- 
nique for comparing stellar populations in clusters that is 
well suited to studies of both cluster and stellar evolution. 



in addition to the interplay thereof. Using a sample of 56 
GCs taken from Piotto et al.'s 2002 HST database, we find 
a sub-linear scaling for RGB stars at the 3-a level, whereas 
the relation is linear for MSTO stars. While the preferen- 
tial self-enrichment of massive GCs, two-body relaxation, 
and evolved BSs could all be contributing to the observed 
sub-linear dependence, further studies with an emphasis on 
selection effects are needed in order to better constrain the 
source of this curious observational result. 
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